Norges Teknisk-naturvitenskapelige Universitet Efficient Time-symmetric Simulation of Torqued Rigid Bodies Using Jacobi Elliptic Functions

نویسندگان

  • Elena Celledoni
  • Niklas Säfström
چکیده

If the three moments of inertia are different from each other, the solution to the free rigid body (FRB) equations of motion is given in terms of Jacobi elliptic functions. Using the Arithmetic-Geometric mean algorithm, [1], these functions can be calculated efficiently and accurately. The overall approach yields a faster and more accurate numerical solution to the FRB equations compared to standard numerical ODE and symplectic solvers. This approach performs well also for mass asymmetric rigid bodies. In this paper we consider the case of rigid bodies subject to external forces. We consider a strategy similar to the symplectic splitting method proposed in [16]. The method here proposed is time-symmetric. We decompose the vector field of our problem in a FRB problem and another completely integrable vector field. In our experiments we observe that the overall numerical solution benefits greatly from the very accurate solution of the FRB problem. We apply the method to the simulation of artificial satellite attitude dynamics.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Norges Teknisk-naturvitenskapelige Universitet Fakultet for Informasjonsteknologi, Matematikk Og Elektroteknikk Hovedoppgave

.........................................................................................................................I PREFACE........................................................................................................................... II TABLE OF CONTENTS................................................................................................... III

متن کامل

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET A method of parametric solution of convolution equations

A variant of the method of moments is developed for parametric solution of convolution equations of the first kind. Two models — the gamma model and the shifted gamma model — are studied in details.

متن کامل

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Parameter Estimation in High Dimensional Gaussian Distributions

In order to compute the log-likelihood for high dimensional Gaussian models, it is necessary to compute the determinant of the large, sparse, symmetric positive definite precision matrix. Traditional methods for evaluating the log-likelihood, which are typically based on Choleksy factorisations, are not feasible for very large models due to the massive memory requirements. We present a novel ap...

متن کامل

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Implementing Approximate Bayesian Inference using Integrated Nested Laplace Approximation: a manual for the inla program

This manual describes the inla program, a new instrument which allows the user to easily perform approximate Bayesian inference using integrated nested Laplace approximation (INLA). We describe the set of models which can be solved by the inla program and provide a series of worked out examples illustrating its usage in details. Appendix A contains a reference manual for the inla program. This ...

متن کامل

Norges Teknisk-naturvitenskapelige Universitet Control Variates for the Metropolis-hastings Algorithm Control Variates for the Metropolis-hastings Algorithm

We propose new control variates for variance reduction in the Metropolis–Hastings algorithm. We use variates that are functions of both the current state of the Markov chain and the proposed new state. This enable us to specify control variates which have known mean values for general target and proposal distributions. We develop the ideas for both the standard Metropolis–Hastings algorithm and...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2006